Quantum Conductance of Achiral Graphene Ribbons and Carbon Nanotubes 
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Explicit expressions of the band spectrum near the neutrality point are derived for armchair and 
zigzag graphene ribbons and carbon tubes. Several spectral features, which were previously observed 
only in numerical calculations, are given an adequate analytic description in terms of elementary 
functions. The obtained dispersion relations are used for a comparison of conductance ladders of 
graphene-based wires; these relations are also beneficial for many other applications. 



INTRODUCTION 

For electrons and holes, graphene ribbons and car- 
bon tubes are one-dimensional wires made of one-atom- 
thick material. In comparison with the two-dimensional 
electron gas counterparts in semiconductor heterostruc- 
tures, the transport of charge carriers in graphene [l[ 
(in particular, quantum conductance [2j) demonstrates a 
number of unusual properties. This paper gives a pre- 
cise analytic description of the conductance of four basic 
graphene wires and specific features of each member of 
the wire family represented in Fig. 1, achiral graphene 
ribbons and carbon tubes. Formally, this problem can 
be considered already "solved" by finding two equations 
that describe the spectrum of a graphene sheet with two 
armchair- and two zigzag-shaped edges [H (in the center 
of Fig. 1Y But spectrum peculiarities near the Fermi en- 
ergy 0, y|, 3, H, [E 0, @1 are far from being obvious from 
the general equations. Here, we show that (i) the band 
spectrum of a metallic armchair (zigzag) ribbon (tube) 
is not the same as for zigzag (armchair) ribbon (tube); 
this difference is given an accurate quantitative descrip- 
tion; (ii) in moving away from the Fermi energy, the bot- 
toms (tops) of conduction (valence) bands in zigzag (arm- 
chair) ribbons (tubes) are shifted towards larger wave 
vectors; and (iii) there exist three types of spectra (con- 
ductance ladders) with equal, irregular, and alternating 
band spacing (ladder step width). By expressing each of 
these features in elementary functions, the understanding 
achieved in previous studies is considerably improved. 



BAND STRUCTURE 

Figure 1 illustrates the parent honeycomb NxAf lat- 
tice and its daughter wire-like structures, armchair and 
zigzag graphene ribbons and carbon tubes, henceforth 
abbreviated as GR and CT, respectively. The lattice la- 
bel indicates that in the armchair direction, graphene 
contains iV hexagons in polyparaphenylene-like chains, 
whereas in the zigzag direction, it has J\f hexagons form- 
ing polyacene-like chains. Hydrogen atoms along edges 
are not shown and not taken into account in the nearest- 
neighbor tight-binding Hamiltonian j^, [l(| ■ 



FIG. 1: Parent graphene NxAT lattice (center) and its daugh- 
ter structures, from lower-left clockwise: armchair ribbon 
(N = oo); (Af,Q) tube (N = oo, 1-st = M + 1-th polypara- 
phenylene chain); (N,N) tube (AT = oo, 1-st = N + 1-th 
polyacene chain); and zigzag ribbon (J\f = oo). In all cases, 
same k-coordinates are used. The cross at (0, |f ) and the 
circle at (0, |j) indicate zero-energy points: x for armchair 
metallic GR, x and o for zigzag metallic CT; the filled circle 
at (771^, §J) indicates the second special point for the arm- 
chair CT. For the zigzag GR, zero-energy point cannot be 
shown on real k x k y plane. 



The 7r-electron spectrum of this model is given by 
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where the hopping integral is used as an energy unit, 
a is the minimal translation distance of the honeycomb 
lattice, and one of the two dimensions, N or TV, is 
supposed to be infinite, implying the continuity of k x 
(for armchair GR and zigzag CT) or k y (for zigzag GR 
and armchair CT). The complementary discrete quan- 
tum numbers are respectively determined by open ends 
and periodic boundary conditions for ribbons and tubes. 
Thus, for graphene ribbons, the k space is (0— 7r,0— tt); 
for armchair and zigzag carbon tubes, the required ex- 
tensions of this space are respectively (0— 2ir, 0— tt) and 

(0-7T,0-27T). 

The spectra of armchair GR (aGR), armchair CT 
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(aCT), and zigzag CT (zCT) are completely determined 
by Eq. ([1]) and by the boundary conditions dictating val- 
ues of the discrete quantum number. An additional equa- 
tion, Eq. comes into play in the discussion of zigzag 
graphcnc ribbons. 

We consider related GR and CT pairs separately. The 
focus is on the energies not far away from the point of 
neutrality, which coincides with zero energy, the Fermi 
energy. In this energy region, the spectrum is described 
by the minus branch of Eq. ([1]) (the minus sign in the 
radicand). Because of the spectrum symmetry, we refer 
only to the conduction one-dimensional bands, that is, 
to the E~ branch with the plus sign in front of the root. 
The valence bands, having the same transverse quantum 
numbers, are just the mirror reflection of the conduction 
bands in the E — plane. The analysis is performed 
for large N and AT. This simplification can easily be 
avoided, but even for N,Af > 10, it is sufficiently good 
for reasonable estimates. 

Armchair ribbons and zigzag tubes. For both 
types of structures, aGR and zCT, k x is a contin- 
uous variable, < \^3ak x < tt. The wave-vector 
transverse component takes discrete values k y = £j/a, 
£i = nj/(AT+l), j = 1,2, ...,7V and = 2nj/M, 
j = 0, 1, ...,Af— 1 for ribbons and tubes, respectively. 

It is known [and also follows from Eq. (TT])] that an aGR 
(zCT) is metallic if j* = 2{N + l)/3 (j* = Af/3) is an 
integer. In this case, EJ = -,{k x =Q) = 0. Otherwise, aGR 
(zCT) is a semiconductor. The twofold band degeneracy 
of electron states in metallic zCT occurs because there 
are two zero-energy points, (0, and (0, 2£j*); equiv- 
alently, (k x = 0,ak y = 2ir/3) and (k x = 0,ak y = Att/3). 
In Fig. 1, these points are marked by a cross and an open 
circle. For a semiconducting aGR (zCT), the index j* of 
the lowest conduction band can be equal to (2JV + l)/3 
or (2A/" + 3)/3 [(Af-l)/3 or (JV + l)/3]. 

The next conclusion that follows from the analysis of 
Eq. (TT]) is that near the cross point for aGR and the cross 
and open-circle points for zCT, the band spectrum can 
be represented as 
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where A zGT = 2A aGR = v^r/Af, X = ^a^/A aGR(zCT) , 
and \i = 0, 1, ... << N '. The validity of these dispersion 
relations is ensured by the condition Afak x << /i =/= 0. 

The upper and lower rows in Eq. ([2]) respectively refer 
to the metallic and semiconducting aGR (zCT). Disre- 
garding the term linear in /i in the upper row yields the 
previously suggested expression for the band spectrum 
of (A/ - , 0) zigzag carbon tubes [B| and armchair graphene 



ribbons [6]. However, as can be easily checked by direct 
calculations of dispersion curves from exact equation {TJ) , 
this term considerably improves the quality of the ap- 
proximate description. 

Armchair tubes and zigzag ribbons. The spec- 
tra of these structures are to be considered separately 
because the band spectrum of the zigzag GR (zGR) is 
determined by two equations. An analytic description of 
the zGR spectrum was repeatedly attempted before but 
never succeeded. 

Armchair tubes. For armchair tubes, ak x = K v j\[?> 
plays the role of the transverse quantum number, k v = 
2m/ /N, v — 0, 1, N—l, whereas < ak y < tt is a con- 
tinuous variable [3j. Equation (fT]) can then be rewritten 
as 
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With an exception of the v — band, all other bands 
arc twofold degenerate. Close to the zero energy, that 
is, for \q\ << 1, q = ak y — 27r/3, and ttv/N « 1 or 
ir(N~v)/N « 1, Eq. is well approximated by 



E-{q) = 
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showing that the spectrum consists of a set of pseudo- 
parabolic bands with bottoms 
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^A aCT « 1, (5) 



^i^; here, A aCT = ir/N. The latter 
quantity would be the band spacing if the term linear 
in q in the radicand in were disregarded, e.g., as in 
Ref. Q. Also worth noting is that for J\f = \/3N, A aCT 
= A zCT = 2A aGR . 

Hence, the spectrum of an aCT contains a nondegen- 
erate band v = with a linear dispersion, followed by 
a manifold of degenerate v and N — v bands with a 
pseudo-parabolic dispersion. This is similar to the case 
of a metallic zCT. But an important distinction is that 
in an aCT, there are two propagating states which have 
different wave vectors, q < qf; CT and q > qf; CT (both 
values of k y are positive) but correspond to the same en- 
ergy. Another difference from a metallic zCT is that the 
obtained correction to the band spacing, 



E^-Er = A a 
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is not linear but quadratic in v/N. This correction is 
small and can therefore be disregarded. 

Zigzag ribbons. Cutting an armchair carbon tube along 
the zigzag direction and "healing" damages of three- 
coordinated cp 2 bonding by hydrogen atoms gives a 
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FIG. 2: One-dimensional bands v = 0, 1, 7 of an (N, N) CT 
(dashed) and v = 0, 1, 14 of an iV-wide zigzag GR, iV = 100 
(solid lines). Curves obtained from exact calculations for aCT 
and zGR are paired with their (often indistinguishable) ap- 
proximations indicated by dotted lines. The linear dispersion 
and all curves in the middle of shaded areas demonstrate the 
aCT band structure Q. The dotted curves correspond to 
Eqs. 0, © and (|10|) . The band minima calculated in accor- 
dance with Eqs. ([5} and (lll[l are respectively shown by circles 
and triangles. 



zigzag graphene ribbon. An essential feature of zigzag 
ribbons is that discrete values of < ak x — n J / \/3 < tt 
are Aiy-dependent. These values can be found by solving 



the equation Q 
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where the minus (plus) sign in the right-hand side corre- 
sponds to the plus (minus) branch in dispersion relation 
(P). By combining Eqs. |T]) and ©, 
As a result, we obtain 



k y can be eliminated. 
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This equation, where only a half of the spectrum with 
positive energies is presented, substantially simplifies 
finding the band structure of a zGR. In particular, the 
extrema of E~ , as a function of kT , are given by solutions 
of the equation 



N 
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For any value of k y , the minus branch of Eq. ((6]) has N 
solutions. One of these solutions is imaginary, Kq = iS, 
if k v falls into the interval 27r/3 + (^/SN)^ 1 < ak y < tt. 
The energies of such states are within a narrow interval, 
Eq(q) < (2N + 1) _1 , and correspond to states localized 
near the zigzag edges. In this case, and under the restric- 
tion N5 » 1, exact equation ([7]) is well approximated 
by 



Eo(q) 
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The inequality NS >> 1 imposes severe restrictions 
on the allowed magnitude of k y . However, the above 
equation describes the dispersion of edge states within 
a half of the actual interval, tt /6 < q < tt /3. It can 
be shown that for small values of q, [V3(N + 1/2)] _1 = 
q c < q < 2q c , the edge-state dispersion is governed by 
Eq = \/3qcxp (— \/3Nq). For the rest of the interval 
27r/3 - tt, i.e., for 2q c < q < tt/6, exact equations ((6|) and 
([7J) must be used. Differences between the exact solution 
of the problem and the approximate description based on 
the Dirac equation Q are discussed in Ref . Q . 

For q < q c , the dispersion of the lowest energy band in 
the zGR spectrum is similar to that of the v = band in 
an aCT (see Fig. 2). To obtain an analytic expression for 
this part and for higher bands, we must use real solutions 
of Eq. ©, K~ fa (v+1/2)tt/N, p = 0,1, ... , v « N. A 
rather cumbersome procedure of finding these solutions 
and exploiting them in Eq. (|6|) yields 
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As explained above, the case v = and q > is special. 

From Eq. (I10|l . the bottoms of the v > bands are 
determined by 



El- = 
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(ii) 

zGR = 7r 2 ( t .+l/2) 2 
4 V / 3AT 2 



These energies are attained at q = q" v 

v = 1,2 We note that the interband separation in 

zigzag GRs is well approximated by A zGR = A aCT /2 = 
7r(2A r )~ 1 . Hence, Eq. (fTTj) can be replaced by E^ — 
{v + l/2)A zGR . The ratio between A zGR and A aCT is 
exactly the same as for zCT-aGR pairs. A comparison 
of the approximations in ([5]), and PU)) with the exact 
calculations (see Fig. 2) proves that these approximations 
are highly accurate. 

By disregarding the term linear in q in Eqs. ([4]) and 
(fT0|) and by passing to the new energy scale ^\t\ (t 
is the hopping integral), which is more convenient near 
the Fermi energy, we can express the band structure of 
graphene ribbons and carbon tubes in a single line 
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FIG. 3: Conductance ladders for zigzag metallic (left), semiconducting (middle) and armchair (right) carbon tubes with different 



circumferences; energy scale (from bottom to up): 
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on the left panel and similarly on the middle and 



right panels. Markers indicating band opening (band bottom for electrons and band top for holes) are calculated from exact 
dispersion relation (1). For each ladder, the straight line G(E)/Go — g v + \E\ represents an approximate dependence on energy 
in the limit as M — > oo (zCT), N —> oo (aCT). This provides a visual estimate of the accuracy of relativistic approximation 
(12), showing that graphene wires do have irregular (left), alternating (middle), and equidistant (right) spectra. Conductance 
ladders for graphene ribbons have practically the same appearance, except that the step height and width are two times smaller. 



where k is the dimensionless wave vector along the ribbon 
(tube) and the index a specifies the structure and the 
corresponding expression for m™. Specifically, 
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where " m" and " s" extensions in labeling indicate metal- 
lic and semiconducting aGR (zCT), and v = 0, ±1, ... for 
all structures except zGR, where v = 0, 1, ... and the case 
v = 0, k > 27r/3 is exceptional. 



Equation (fT2"|) has the form of a one-dimensional rel- 
ativistic energy-momentum relation in its conventional 
representation with the speed of light equal to unity. This 
is not more than a formal analogy. 

As follows from Eqs. (fl"2"|) and (fl~3"|) . the spectra of 
achiral graphene ribbons and carbon tubes can be clas- 
sified into three groups: (i) metallic spectra with equally 
spaced bands, as for the aCT and zGR; (ii) metallic spec- 
tra with a regularly irregular band spacing, as for the 
zCT and aGR; and (iii) semiconducting spectra with the 
band spacing alternating between m and m (between 
and w)> as f° r semiconducting armchair ribbons 
(zigzag tubes). 

CONDUCTANCE 



In the framework of the Landauer approach [111 . Il2j , 
the zero-bias, zero-temperature conductance of an ideal 
wire is equal to 



(14) 



where Go = 2e 2 /h is the conductance quantum, <?„ is the 
degeneracy of band states (spin degeneracy 2 is included 
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into Go), and the transmission coefficient T„ is zero or 
unity, depending on whether the i^th band is open for 
charge carriers or not. 

For the band structure specified in Eq. (TT2"|) , T V {E) = 
®{E-m%) for conduction bands and T V (E) = ®{\E - 
m"|) for valence bands, where Q(x) is the Heaviside step 
function. Thus, the quantum conductance as a function 
of energy has the form of a ladder, symmetrically as- 
cending with an increase in E for electrons, and with a 
decrease in the energy for holes. The height of the fth 
ladder step is determined by the band-state degeneracy 
g u =o = 1 (2) for zGR (aCT), otherwise, g u = 2 (4) for 
GR (CT). Three types of the band spectrum identified 
above can be translated into three corresponding types 
of conductance ladders with regular (aCT and zGR), ir- 
regular (zCTm and aGRm), and alternating (zCTs and 
aGRs) width of steps. 

Needless to say, the appearance of G(E) illustrated in 
Fig. 3 depends on the energy scale determined by the 
ribbon width (tube circumference). It can vary from the 
ladder-shaped curve, ascending with the increase of \E\ 1 
to a straight line G — Go(g l/ + \E\). The latter depen- 
dence shows that metallic and semiconducting graphene- 
based wires are indistinguishable as classic conductors. 
Furthermore, if the irregularity of the step width is re- 
solved, the difference between conductance ladders of 
metallic zCT (aGR) and aCT (zGR) is not questioned, 
but they look identically in a larger energy scale. This 
can be seen in Fig. 3, which illustrates conductance in 
different energy scales. 

The above consideration concerns four basic graphene 
wires which are discussed in more detail elsewhere 
[l3j |. It provides a useful reference for calculations of 
electron and hole coherent transmission in various types 
of graphene contacts. 
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* This paper is published in Zh. Eksp. Teor. Fiz. 135, 139 
(2009). There is no difference in the presentation of the 
results between the present arX version and published 
paper. However, there are some small but nevertheless 
important textual changes in this version, where a few 
misprints are eliminated. Also, Eq. (7) reads as presented 
here, and not as in the journal version. We regret for this 
inconvenience caused by the editorial process which, in 
our opinion, could have better feadback with the authors 
and better timing. 
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